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Abstract 

The structure of flexible polymers endgrafted in cylindrical pores of diame- 
ter D is studied as a function of chain length N and grafting density a, assum- 
ing good solvent conditions. A phenomenological scaling theory, describing 
the variation of the linear dimensions of the chains with a, is developed and 
tested by Molecular Dynamics simulations of a bead-spring model. 

Different regimes are identified, depending on the ratio of D to the size of 
a free polymer iV 3 / 5 . For D > JV 3//5 a crossover occurs for a = a* = iV -6 / 5 
from the "mushroom" behavior (R gx = Rgy = Rgz = N 3 ^ 5 ) to the behavior of 
a flat brush (R gz = a^ 3 N, R gx = R gy = a" 1 / 12 N 1 / 2 ) , until at a** = {D/Nf 
a crossover to a compressed state of the brush, (R gz = D,R gx = R gy = 
(A /_3 Z)/4cj) 1 / 8 < D), occurs. Here coordinates are chosen so that the y-axis is 
parallel to the tube axis, and the z-direction normal to the wall of the pore 
at the grafting site. 

For D < iV 3 / 5 , the coil structure in the dilute regime is a "cigar" of 
length Rgy = ND~ 2 / 3 along the tube axis. At a* = (7V.D 1 / 3 ) -1 the structure 
crosses over to "compressed cigars" , of size Rgy = (aD)- 1 . While for ultrathin 
cylinders (D < N 1 ^) this regime extends up to the regime where the pore is 
filled densely (a = D/N), for N 1 / 4 < D < N 1 / 2 a further crossover occurs 
at a*** = D- g / 7 N-^ 7 to a semidilute regime where R gy = (N 3 D/Aa) 1/8 still 
exceeds D. For moderately wide tubes (N 1 / 2 < D < iV 3 / 5 ) a further crossover 
occurs at a**** = N S D~ 7 , where all chain linear dimensions are equal, to the 
regime of compressed brush. 

These predictions are compared to the computer simulations. From the 
latter, also extensive results on monomer density and free chain end distribu- 
tions are obtained, and a discussion of pertinent theories is given. In particu- 
lar, it is shown that for large D the brush height is an increasing function of 
D-\ 



1 Introduction 



Polymer brushes are layers of flexible linear polymers on substrate surfaces produced 
by endgrafting one chain end via a special chemical endgroup of the macromolecule 



of colloid stabilization, lubrication, wetting, adhesion, etc. and hence they have 
found abiding interest. Polymer brushes also pose challenging theoretical problems, 
since the confinement due to the grafting at the substrates causes the emergence 
of many distinct characteristic lengths needed even to describe the structure of a 
single chain in a brush [01[3IBj, and although the density profile of the monomers 
in the direction normal to the grafting surface has been studied since a long time 
ago P HU1 HD 112 13 13 113 CHI El, the degree to which this profile is understood 
quantitatively still is a subject of current discussion [TH] . 

Note that the overwhelming part of the very rich literature on polymer brushes 
(see [U 121 El IH E] for reviews) deals with brushes on planar flat surfaces only. Occa- 
sionally the problem of brushes on the outside of spherical particles {e.g. [3I2D1P1J} 
and on the inside of spherical cavities j20J 1221 l2Sj wa s considered (as well as weakly 
curved flexible membranes, e.g. fM\ I25j). 

However, relatively little work has been done on grafted flexible polymers inside 
cylindrical pores (apart from a few comments in ffiH 123] and a study of a single 
grafted chain in a tube [2H], a problem that is closely related to the behavior of 
single non-grafted chains confined in tubes [23 12H| ) ■ We note that small mushrooms 
grafted to cylindrical pores were considered as a means to control electroosmotic 

flow |22| • 

In the present paper we intend to contribute to fill this gap, by presenting a 
study of polymer brushes grafted on the inside of cylindrical pores of radius R, 
varying also the chain length N and grafting density o over a wide range. Such 
systems should be possibly also relevant in a biophysical context (cylindrical pores 
in biological membranes, blood vessels or other transport pipes or tube-like objects 
in biosystems, in which proteins or other biopolymers can be adsorbed), but this is 
outside of consideration here. Rather this paper is devoted to an understanding of 
the characteristic lengths that describe the configurations of the endgrafted chains 
in such a cylindrical confinement. In Sec. II, we formulate a (somewhat qualitative) 
scaling description, pointing out also the relation to the problem of non-adsorbed 
chains in cylindrical pores at various volume fractions of the monomers [23 EDI EH1 
E21 EH! El| • In Sec. Ill, we describe the model used for our simulations to test this 
theory, while Sec. IV presents our numerical results for the linear dimensions of the 
chains. Sec. V describes results for the monomer density distribution across the 
channel, as well as corresponding results for the distribution function for the free 
chain end. Sec. VI then summarizes our conclusions. 

2 Theoretical background: Scaling considerations 

A very brief discussion of polymer brushes in cylindrical geometry was given in 
Appendix C of [2B], however, many of the characteristic lengths relevant for our 
study were not given, and also a detailed derivation of the diagram describing the 
different regimes in the plane of variables R and a was not presented; we give such 
a justification here, in order to make clear which underlying assumptions are made, 
and in order to be able to derive the new results on various characteristic lengths 
that are of interest here. 

2.1 Mushrooms and Polymer Brushes on Flat Substrates 

To set the scene, we recall the "blob model" of Alexander |J] and de Gennes [TU] 
and its extension to deal with a nonuniform monomer concentration within a brush 
[0] , see Fig. ^ For simplicity, all lengths will be measured in units of the effective 
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N segments simply is \fN , and both grafting density a and volume fraction <fi are 
treated as dimensionless quantities. Only the good solvent regime is treated, and 
- as usual in the spirit of scaling approaches [35] - prefactors of order unity are 
suppressed throughout. Also the exponent v describing the linear dimensions of a 
self-avoiding walk in d = 3 dimensions is simply approximated at the Flory value 
36J v = 3/5, rather taking the more accurate value v ~ 0.588 [SZ1 EH!- Thus an 
isolated polymer chain endgrafted at a planar wall, a "polymer mushroom" , has the 
same gyration radii components as in the bulk, 

Rgx = Rgy = Rgz = N 3/5 . (1) 

Here, and in the following, we orient the z-axis normal to the grafting surface while 
the x and y-axes lie in the grafting plane (in the case of a a planar brush). 

Considering then grafting densities a large enough such that the individual mush- 
rooms would significantly overlap, we recognize that (for a <C 1) we obtain a semidi- 
lute polymer brush, in which the screening length £ is of the order of <r~ 1//2 , which 
then is taken as the blob-radius in Fig. ^ It is implied, that inside a blob excluded 
volume statistics persists, only n monomers from one chain are contained inside of 
such a blob, and no monomers of any other chains. Since a -1 / 2 = n 3//5 , in analogy 
to Eq. we conclude that a blob contains n = <t _5//6 monomers, and every chain 
contains iV Wo6 = N/n = a 5 ^ 6 N blobs. According to the Alexander j^j-de Gennes jTU| 
picture, the chains in the brush form "cigars", i.e. they form linear chains of blobs 
stretched out along the z-axis (fig. HI upper part), and hence the brush height h is 
of the same order as R gz , namely 

h = R gz = £N blob = cr^N. (2) 

Taken literally, this "cigar" configuration would imply that R gx = R gy = £ = cr^ 1 ^ 2 , 
which in reality is not the case: rather the chains in x and y-direction have a structure 
resulting from a random walk of blobs, i.e. 

Rgx = Rgy = = a~^N^. (3) 

While the Alexander [0]-de Gennes jTUj picture also implies that the monomer den- 
sity profile 4>(z) is essentially uniform for z < h, and the free chain ends are all in 
the outermost blob, i.e. in the region h — £ < z < h, the structure in reality is 
much more complicated [13 El El EH ECU EH ESI EH E3 EEJ • Rather than staying con- 
stant, p(z) smoothly decreases with increasing h (and in the strong stretching limit 
approximately is described by a parabolic variation, <j>(z) = 0(0) [1 — z 2 /h 2 }). Due 
to this decrease of the monomer fraction <j>{z) with the distance z from the grafting 
surface, also the screening length £ is not a constant, but increases according to 
£(z) = [^(z)] -3 / 4 [33]. From a self consistency consideration one can then estimate 
the diameter dfi of the final blob as 6j dfi — N 2 ^ . All these concepts have been 
verified by corresponding computer simulations [UJ El EB EUl EH EH] • 

2.2 Cylindrical Tubes 

We now consider a cylindrical tube of diameter D and length L. We orient the tube 
axis in y-direction and choose in this direction periodic boundary conditions. 

We first note that there is a simple geometric relation between the grafting 
density a of anchor points at the cylinder surface and the average volume fraction 
of monomers contained in the cylinder. Since the area of the cylinder surface is 
A = ttDL, the total number of monomers M contained in the cylinder is 



The cylinder volume being V = ttLD 2 /4:, we find 



cf) = M/V = ANa/D (5) 

The condition that the maximum volume fraction that can be reached in a cylinder 
corresponds to a dense melt, i.e. = 1, implies the following constraint for physically 
meaningful grafting densities a 

a < D/(AN). (6) 

We first discuss the regime of extremely low grafting densities. For wide tubes, 
D > iV 3//5 , we still have mushrooms, i.e. Eq. still holds. A more interesting 
regime occurs for narrow tubes, 

1 < D < iV 3/5 , (7) 

where the chains form "cigars" stretched along the y-axis (Fig.|2J). The situation for 
a — > closely resembles the case of non-grafted chains confined to very narrow tubes 
under good solvent conditions [271 |2H]- One has R gx = R gz = D, while the longi- 
tudinal ^/-component R gy is much larger. Assuming that the "cigar" configuration 
can be treated as an array of iVf,^ = N/n blobs of diameter D, with n monomers 
per blob, D = n 3//5 , we readily obtain 

R gy = DN blob = DN/n = ND~ 2 ' 3 (8) 

As it should be, for D = N 3 ^ 5 this yields a smooth crossover to R 9y = iV 3//5 . 

Next we estimate how small the grafting density a has to be in order that Eqs. (JIJ 
or (jHJ) hold, respectively. The corresponding boundary in the D — a plane is obtained 
by the condition that the grafting density is such that the mushrooms or cigars just 
touch each other. 

In the case of mushrooms, D > iV 3 / 5 , this boundary is reached for an overlap 
grafting density a* given by 

a* = N~ 6/5 , D > N 3/5 (9) 

since each mushroom covers an area of order N 6 ^ 5 when one projects the monomer 
coordinates of the chain onto the grafting surface. 

In the case of cigars, elongated along the y-axis, D < iV 3//5 , we simply have to 
equate the density in the tube to the density inside a cigar, using Eq. © 

0* = N/ (D 2 R gy ) = D~ 4 ^ 3 . (10) 

This overlap density 0* is easily converted into an overlap grafting density a* using 
Eq. ©, 

a* = D(j)*/(AN) = l/(4iV£> 1/3 ), 1 < D < N 3/5 . (11) 

Also in this case we note that for D = N 3 ^ 5 both expressions Eqs. Q, (fTTj) for the 
overlap grafting density are of the same order, as it should be. 

We now consider the behavior for a > a*. In the regime where D > N 3 ^ 5 a 
crossover from mushroom to brush takes place, similar to the case of planar grafting 
surface. To leading order, Eq. (j2J) still holds. However, a strong compression of the 
brush sets in when the brush height h (and thus also R gz ) become comparable to 
the tube diameter. This yields a crossover grafting density a** separating a weakly 
compressed brush (described still by Eq. Q) from a strongly compressed brush, 
with R gz = h = D, for a > a**, 



In the regime of the compressed brush, the density of monomers is approximately 
uniform throughout the tube. Then the lateral linear dimensions are, similar as for 
chains in a semidilute solution in the bulk (£ = 0~ 3 / 4 ), given as 

R 9X = R gy = £JVJ2 = 0~ 1/8 iV 1 / 2 = (4a)- 1 / 8 J D 1 /s iV 3/8 ) a > a ** (13) 

At the boundary, a = a**, Eq. (fTHj) yields linear dimensions of order 

R gx = R gy = D-^N 3 ^ ,a = a*\ (14) 

and hence there occurs then again a smooth crossover to Eq. Q. 

Finally, we address the behavior for o > a* in the narrow tube regime D < iV 3 / 5 . 
While for a = a* we have found "cigars" (Fig. |2J) as described by Eq. (jHJ) that just 
touch at their ends, we now expect "compressed cigars" in which the excluded 
volume interactions are partially screened. The polymer chains overlap themselves 
and near the cigar ends, neighboring cigars overlap each other. The screening length 
(blob radius) £ then is smaller than the tube diameter, unlike the situation depicted 
in Fig. |21 namely 

£ = ^ = (4Na/D)-V\ (15) 

where in the last step Eq. (jHJ was used. 

For the range of grafting densities a* < a < a*** (a*** will be estimated below) 
we have compressed cigars with linear dimensions R gx = R gz = D and linear dimen- 
sions R gy > D along the y-axis. Assuming that the volume D 2 R gy is of the same 
order as the total volume taken by all the blobs of a chain, ^ 3 Nu b = ^£ 4//3 , one 
finds 

Rgy = ZT 2 iV£ 4/3 = D^Ncj)- 1 = l/(4aD) (16) 

The result R gy = D^Ncf)^ 1 is identical to the corresponding result for non-grafted 
chains confined to narrow tubes [HOI El]; as expected. In this quasi-one- dimensional 
situation the "cigars" are laterally compressed, each chain overlaps itself to some 
extent (excluded volume is respected only up to the length £, Eq. (|15|)). but neigh- 
boring chains are still essentially segregated from each other and overlap only near 
the ends of the "cigar" . 

Now the boundary a*** of this regime can simply be estimated from the condition 
that Rgy, Eq. flTS]) , equals the result in the compressed brush regime, Eq. (|T3)). which 
yields 

4 

Note that for D = N 3/5 this boundary starts at a*** = (l/4)iV~ 6/5 , i.e. it splits off 
from the boundary a* in the (D, a) plane (Fig. EJ). The other end of this boundary 
is reached for the condition 0=1, i.e. 4cr*** = D/N (Eq. (JSJ)) and hence D = N 1 ^. 
This means, for ultra-narrow tubes (defined by the condition D < N l l A ) Eq. (fT6*|) 
holds from a = a* up to a = D/4N, so the boundary a*** cannot be reached. 

For tubes that are not ultra-narrow (D > N 1 ^) one encounters at a*** a crossover 
to a regime of semi-diluted cigars where now two of the three linear dimensions 
R gx = R gz = D, while the third linear dimension, R gy , is given by Eq. (fT3|) . However, 
the cigars now start to overlap more and more, until at a grafting density a**** the 
"natural size" R gy of a chain in a semidilute solution of density <fi reaches D, 

R gy = ^/j; = 0- 1 / 8 iV 1 /2 j a *** <a< (j**** . (18) 

While in the regime described by Eq. (|T5|) R gy still exceeds D, and therefore the 
linear dimension of a chain in y-direction exceeds the other linear dimensions, we 



denote this regime as "overlapping cigars" . The grafting density a**** is found from 
the condition that R gy = D, and hence 

D = £Nlll = (ANa/D)- 1/8 N 1/2 , a**** = ^N 3 D~ 7 (19) 

For a**** < a < D /AN we enter the same region of compressed brush that we have 
already encountered for D > N 3 ^ 5 for a** < a < D /AN. The line a**** starts in the 
(D,a) plane also for D = N 3 ^ 5 at the point a = jN~ 6 ^ 5 , i.e. at a*, and ends for 
(j) = 1(<7 = D/AN) at the diameter D = N 1/2 (Fig. EJ). 

The behavior of the brush in a tube becomes somewhat more transparent if we 
discuss it in terms of the scaled variables D = D/N 3 ^,a = o~iV 6//5 , since then (in 
the dilute and semidilute regime) the dependence on the chain length is absorbed 
completely. Fig.|3]shows qualitatively the dependence of the chain linear dimensions: 
For D > 1 the scaled radius R gz = R g /N 3 ^ crosses over from unity at a = 1 to the 
power law cx 1//3 and at a = a** = D 3 to D, while the other components of the radius, 
R gx = Rg X /N 3 / 5 , Rgy = R gy /N 3 ^ cross over from unity at a = 1 to the power law 
cr -1 / 2 and at a** (where R gy = Z) -1 / 4 ) a further crossover to (D/a) 1 ^ 8 occurs. 

Similarly, for D < 1 the radius R gy is D~ 2 / 3 for small a and crosses over at 
a* = D^ 1 / 3 to a power law (cx-D) -1 . For D < iV~ 4//15 this is the only crossover that 
occurs, while for larger D additional crossovers are found at a*** (and a****, if D 
exceeds N~ l l w ) with a*** = D' 2 , a**** = D' 7 . 

For a > a*** we have R gy = D, but only at a — a**** the radius of a chain in a 
bulk semidilute solution of density <fi corresponding to a would be of order D and one 
goes over to compressed brush behavior where R gy decreases to linear dimensions 
smaller than D, namely R gy = (D/a) 1 ^ 8 , as above. 

Note that in the case D < 1 the behavior of the brush in the tube in many aspects 
is similar to the behavior of non-grafted chains in a tube, as discussed by (HOI EH] • 
One can see this converting a to 0, using Eq. (J3J); however, for the latter problem 
the crossover at a**** has a different meaning: then all three linear dimensions 
Rgx, Rgy, Rgz are smaller than D, given by R ga = (D/a) 1 ^ 8 , for a = x, y,z while 
for a < a**** we have R gx = R gz = D for the non-grafted confined chains. For 
the grafted confined chains, we also have R gx = R gz = D in this regime, but R gx 
and R gz may differ from each other by factors of order unity (which are suppressed 
here in our scaling description), while for the non-grafted chains these components 
are identical by the symmetry of the problem. For a > ^****, however, there is a 
difference on a scaling level, since the grafting of the chain ends at the walls creates 
an elastic energy in the brushes, causing chain stretching in the normal direction 
(z), and R gz stays at its (maximal) value which is of order D. Finally, we note that 
in the scaling diagram of Manghi et al. the line corresponding to a*** was not 
included. From Fig. HI however, it is evident that an intermediate regime between 
the isolated cigars (R gy = D~ 2 ^ 3 ) and the semidilute regime (R gy = (D/a) 1 ^ 8 ) is 
mandatory to provide a smooth crossover of the linear dimension R gy and a*. Note 
also, that on the crossover line a**** where R gy = D and hence all linear dimensions 
are of order D, the behavior of R gx changes: we have R gx = D for a < a**** but 

Rgx = Rgy = (D / d) for O > 

3 Model and some Remarks on Simulation Meth- 
ods 

Being interested in phenomena of mesoscopic length scales of order lOnm (or more) 
where chemical detail only enters into prefactors of various quantities, and a scaling 



where atomistic details are not included is mandatory for the sake of computational 
efficiency of our computer simulation |39] I3HJ El G21 EH- There is compelling evi- 
dence that the generic properties of flexible polymers in solution and melt can be 
described by simple bead-spring models where each bead may represent a few chem- 
ical monomers along the chain backbone |39] 147] . and actually such models have been 
used for previous simulations of polymer brushes [SI EJ E3 CHI HH1 EH] successfully. 

In this bead-spring model the interactions between the beads modelling the 
monomer units are described by a Lennard- Jones (LJ) potential that is truncated 



at r c and shifted there to zero, 



U LJ {r) = As 



(j \ 12 /rr \6 



127 

+ £ S096- r<r = (20) 



where r c = 2 7 ^a mm . Here e characterizes the strength and a mm the range of the 
LJ-potential. Note that the minimum of Ulj occurs at r m j n = 2 1 ' 6 a mm , without the 
additive constant £^§q the depth of this minimum would be Ulj{t = r min ) = —e. 
Due to the shift of Ulj at r c = 2r min , the potential in Eq. (|2U|l is continuous at all 
r > 0. 

In addition, monomers bonded to each other as nearest neighbors along a polymer 
chain interact with a finitely extensible nonlinear elastic (FENE) potential |33] 

R n \\ A r 2 



Ufene{t) = -15e ( log I 1 - — f ) , R = 1.5a mm . (21) 



The choice of these parameters ensures that the minimum of the total potential 
between two bonded monomers along the chain occurs for [H^HO] r ~ 0.96cr mm , dis- 
tinct from r m i n w 1.12<r mm . This misfit between the two distances ensures that there 
is no tendency of monomers to form a simple crystal structure even if the density 
is very high and/or the temperature is rather low. However, in the present paper 
we shall consider only temperatures in the good solvent regime of our model. This 
means, the temperature T must exceed the B-temperature, which for the present 
model has been established as (17] ksT je « 3.3. Thus, we choose ksT/e = 4 for our 
simulation. To simplify the notation we choose ks = 1 and take e = 1 and a mm = 1 
as our units of energy and length, respectively. 

The grafting wall is represented by atoms forming a triangular lattice, wrapped 
around a torus by a periodic boundary condition. This choice of wall is inspired 
by carbon nanotubes. However, the potential with which effective monomers and 
wall atoms interact is chosen to be the same as the monomer-monomer potential, 
Eq. (|20p. Also, the wall atoms are fixed in the rigid positions of an ideal crystal 
lattice. This guarantees that no monomer can cross the wall of the tube and get 
out of it. Of course, the coarse-grained nature of our polymer model would make a 
chemically realistic description of a cylindric pore meaningless. The polymers are 
grafted hexagonally on a cylinder with radius R— 1.0, thus defining a grafting density 
a while the tube atoms are placed a little further - on a cylinder with radius R. It 
might be more realistic, of course, to choose randomly distributed grafting sites, 
however, the experience with flat brushes ^7] shows that it makes little difference 
whether the grafting sites are distributed regularly or at random. The size of the 
simulation box along its axis y is L y and the number of static particles, forming the 
tube Ntubei are different for simulation of tubes with different radii R: for R = 15, 
L y = 81.62, N tube = 1024; for R = 21, L y = 76.18, N tube = 2400; and for R = 31: 
L y = 67.47, N tube = 1240. 

Simulations are performed by a standard leapfrog algorithm at constant temper- 
ature, maintained by a Nose- Hoover thermostat [48J, combined with a simple "hot 



"particle temperature" differs from the average value too much. Choosing the mass 
of an effective monomer as m = 1, the MD time unit is to — { (7 mrn m /^ £ ) 1 ■> an d 
the time step of integration 5t = 0.0005to- 

Particular care has to be devoted to the construction of the initial configuration 
of the grafted polymers and their equilibration. The construction of the starting 
configuration for given (N, a, R) begins with the growth of the chains in a tube with 
radius Ri n u — 2-R and N = N max , where N max is the number of monomers in the 
longest chain, investigated for these a and R. After placing the first monomer at 
the grafting site, each consecutive monomer is placed at a distance 0.75_R from the 
preceding one along a direction r = 2.0 * e 0z + e rand where e 0z is a unit vector, 
pointing radially from the grafting site to the tube axis, and e rand is a unit vector 
oriented at random. After an initial equilibration of this system for about 10 5 MD 
steps, we rescale the system radius to the final system with desired radius R and 
again let the system equilibrate for 10 5 steps. From the final configuration of this 
system, by cutting systematically the polymer chains to different values of N, we 
obtain the "starting configurations", which we additionally equilibrate for about 
3.10 5 steps and (depending on iV and a) then collect statistics during the following 
2.10 5 — 4.10 5 integration steps. Typical examples of equilibrated brush configurations 
in this cylindrical geometry are shown in Fig. For all three tube diameters D = 
30,42, and 62 as well as for chain lengths 8 < iV < 64 we then sample brush 
profiles and components of gyration radius, R g , and end-to-end distance vectors of 
the chains, R e , for seven different grafting densities a. In our simulations we have 
refrained, however, from creating elongated cigar-like configuration, cf. Fig.^; since 
the distribution p{y) of R e is expected to be bimodal (that is, the center-of-mass 
of the chains is either to the left, or to the right of the grafting site, as shown in 
Fig.[2b)- The equilibrium p{y) is expected to produce a pronounced minimum above 
the grafting site at y — and, therefore, a well equilibrated distribution p(y) would 
imply that the chain end passes the coordinate y = many times which can hardly 
be possible within the typical time interval of a MD simulation. 

4 Linear dimensions of chains in brushes in cylin- 
drical pores 

In order to test our model and check that for the rather short chain lengths used in 
the present work (typically 8 < iV < 64) results are obtained that one can compare 
to the theoretical predictions, which really address the asymptotic regime N — > oo 
only, we first present data for brushes on flat substrates in Figs. El One recognizes 
that for small a, such as a ~ 0.0115, the linear dimensions in the z-direction, per- 
pendicular to the substrate, and in the ^-direction, are rather similar for the range 
of chain lengths that is studied, but the effective exponents u e ff in the relations 
(Ryy) oc N 2Ue ff , (R 2 gz ) N 2Ve ff somewhat differ from the theoretical value in the 
mushroom regime, v = 3/5 (Eq. For the grafting density o m 0.1039, how- 

ever, we recognize that the scaling behavior of the chain linear dimensions strongly 
depends on direction, (R 2 Z ) oc iV L86 while (R gy ) oc iV L16 , and approaches the the- 
oretical behavior that is expected for much larger N, (R 2 X ) oc iV while {R 2 gz ) oc N 2 
(Fig. Et). Since comparatively short chains are used, it would be premature to ex- 
pect that one reproduces the asymptotic exponents of the scaling theory precisely. 
Also the dependence of both {R 2 gy ) and {R 2 gz ) on a is roughly compatible with the 
predicted power laws (Eqs. ©, ©)• Of course, for polymer brushes on flat sub- 
strates it would be easier to equilibrate significantly longer chains: e.g., jlH] used 
50 < iV < 200 for a closely related model, and ^H] used 16 < N < 512. Even then 



However, for the brushes in the cylindrical pores good equilibration is signifi- 
cantly more difficult; therefore we restrict the present work to relatively short chains. 

Thus, it cannot be a surprise that in the scaling plots (Figs. Ek,b) a perfect 
collapse on master curves is not obtained, and rather some systematic deviations 
due to corrections to scaling are clearly apparent. Such corrections to scaling are 
presumably also responsible for the fact that the effective exponents u e ff in the 
mushroom regime of Figs. \§jp,c deviate from their theoretical value. Note also that 
in the present work a somewhat more restricted range of grafting densities is studied 
than in [IB], where much lower grafting densities were included to ensure to reach 
the mushroom regime fully. Our data in Fig. |U] even for the smallest values of a are 
already affected a little bit by the crossover towards the brush regime, since this 
crossover occurs rather gradually. 

In ^H] also the problem was discussed that for temperatures not too far from the 
Theta temperature the excluded volume interaction on the scale of a blob diameter 
may be so weak that then in the regime of high grafting densities inside of a blob 
gaussian chain statistics still holds. As a consequence, then a further crossover that 
occurs for large a under good solvent conditions has been identified. However, this 
crossover gives rise to rather small and subtle effects, and hence can be safely ignored 
in the present context. 

We now turn to our data for brushes in cylindrical pores. Figs. IHllTOl present 
the linear dimensions (R gz ) and (R 2 gy ) for three choices of the pores diameter, D = 
30, 42, and 62. For the narrow pore (D = 30), one sees a weak increase of (R gz ) 
with a for small N, while for larger N no longer any increase occurs. This behavior 
is in marked contrast with the behavior of the brush on the flat substrate. In 
contrast, (R 2 gy ) decreases with a monotonously, while the theoretically predicted 
double crossover (Fig. 0J) is not seen. For the wider tube, such as D = 62, one still 
sees a more pronounced increase of (R gz ) with er, but if one would fit a power law to 
this increase, the effective exponent would be distinctly smaller than the theoretical 
value 2/3, that should still be observable in an intermediate regime of er, see Fig. HJ 
However, since the crossovers (from the mushroom regime to the brush and from 
the brush to the compressed brush) do not occur at sharply defined values of a, 
but are rather expected to be smeared out over a range of values a that may be a 
full decade or even more. Therefore instead of the two "kinks" in Fig. 0] a smooth 
gradual increase with a is noted. Similarly, when we study {R 2 gz ) as a function of 
N (FigJTUfc), one can still see straight lines on the log-log plot, but the observed 
effective exponents are somewhere in between 2v = 1.2 and 2, indicating the onset 
of compression due to the confinement of the brush in the tube. However, such 
effective exponent which simply reflect slow crossovers do lack any deep physical 
significance, of course. 

Motivated by the scaling analysis of Sec. II, we discuss the extent to which the 
chain length dependence can be absorbed in scaling plots of R gy and R gz versus a 
(Fig. ITT)) . Of course, only for D > 1 the first crossover from mushrooms to the 
almost unperturbated brush is expected to stay unaffected by D. Since varying N 
at fixed D means that different values of D occur, the positions of the various other 
crossovers in Fig. do depend on iV via the iV-dependence of D. So we expect 
that the curves for R 2 Z vs a at large enough a start to split off from the master 
curve that describes both the mushroom and brush regime and saturate at plateaus 
corresponding to the different values of D. Some evidence for the splaying out of the 
curves on the right hand side of the scaling plot for R gz can indeed be seen in Fig. fTTl 
However, for the parameters accessible in our study the crossovers at er* = 1 and 
at a** = D are too close together, so the intermediate regime where R = a 2 / 3 is 
never clearly realized. In fact, in order to clearly separate these crossovers, one would 



small iV corrections to scaling come into play also and obscure the behavior further. 
Part of our data falls into the region D < 1, where a different pattern of crossover 
behavior applies, but again the different regions are clearly separated only if D <C 1 
is achieved. We could not reach this region due to difficulties of equilibration. 

Thus Fig.Eldoes not give strong evidence for the scaling theory of Sec. II: clearly 
these scaling considerations are useful only for the special limit where both N and 
D are very large, a being very small, so that a can be varied over many decades, 
as well as D. This limit is not easy to reach in our simulations. Nevertheless the 
scaling theory gives a useful orientation: note that the data for the scaled quantities 
Rg and R 2 gz in Fig. El show remarkably little variations in spite of the fact that 
both N and a are varied over about an order of magnitude each. 



5 Monomer density and chain distribution across 
a cylindrical pore 

We start this section again with a brief glimpse on the behavior of our model in 
the limit of a flat grafting surface - Fig. El As expected, one sees a smooth decay 
of the monomer density profile <p(z) with increasing distance z from the surface. 
When N increases, the height h of the brush (region where <p(z) is appreciably 
nonzero) strongly increases (Fig. fl2*k). Since from Eq. (J2J) we expect that h oc N, 
it is interesting to plot <j>(z) against a rescaled coordinate z/N (Fig. [T2b). One 
can see that for large N the curves indeed converge smoothly towards a master 
curve, and this master curve resembles the simple parabolic shape predicted by the 
selfconsistent field (SCF) theory in the strong stretching limit ^21 El El 

cj)( z ) = 0(0) [1 - (z/h) 2 } N -> oo, a -> 0, Na 1/3 finite (22) 

Eq. fl2H) implies a nonzero intercept 0(0), while the actual (f>(z — > 0) — > 0, of course 
(Fig. H2J): so the 0(0) of Eq. (}22j) is interpreted as an extrapolation of the flat part 
of the profile; also <fi(z) = for z > h, which for finite N is not true. However, 
numerical implementations of SCF theory that do not invoke the strong stretching 
limit, such as the Scheutj ens- Fleer theory |5T| . can describe the smooth decay of the 
tail of the profile <j>{z) for z > h 52\, but this is not attempted here. In order to 
provide a precise definition of h from the observed profiles <p(z), such as shown in 
Fig. El we follow previous practice (e.g. [TSj) to use the first moment (z), 

oo oo 

&=!<*), (z) = f z<j>(z)/ J <j>{z)dz. (23) 



It is easily checked that Eqs. (|2*2*J) and f!23[) are compatible with each other, but 
Eq. f!23j) is generally applicable. The insert of Fig. El gives evidence that the scaling 
regime h (x N has not yet fully been reached by our rather short chains, but at least 
there is evidence that the chains are stretched away from the surface. Similar data 
have been generated by many previous simulations on various models jUl El El 
El Ej- The delta-function like peak for the first monomer (at z = 1.0) simply is 
due to the fact that the wall atoms forming the flat surface are at a fixed distance 
z = 1.0 away from the first monomer of the chain. 

We also draw attention to the weak oscillations in the density profile ( "layering" ) 
that is visible for small z. Such effects arise from the "packing" of monomers near 
the grafting surface, which resembles a hard wall. Such packing effects are neither 
described by the scaling theory of Sec. II nor the SCF theory, of course. 

Also the distribution function p(z) of the end monomers is of interest (Fig.Efc). 



prediction of the Alexander [Sj - de Gennes picture, according to which the chain 
ends are all localized at a distance z ~ h from the grafting surface. According to 
the strong stretching limit of SCF theory, however, one would predict a distribution 
that is spread over the entire brush but vanishes at z — h with a divergent slope 

<f> e {z) oc z[l - (z/h) 2 ] 1/2 , N -> oo, a -> 0, a 1/3 iV finite (24) 

From Fig. IT2"b it is clear that our data do not reach the strong stretching limit of 
the SCF theory, however. 

For polymers grafted to a flat wall it would be possible to simulate much longer 
chains, of course, and, in fact, in ^Bj for a closely related model (perfectly flat 
wall with no corrugation whatsoever and using an integrated Lennard- Jones-type 
repulsive wall potential) chain lengths up to N = 512 were studied. However, even 
for the largest chains the strong stretching limit of the SCF theory was not reached. 
Since the main interest of the present paper concerns chains in cylindrical tubes, and 
for this case it is probably not feasible experimentally to graft extremely long chains 
at the tube walls, we shall not emphasize the comparison with the strong stretching 
limit of the SCF theory in the present paper; just as the scaling considerations of 
Sec. 2 it is useful for a first orientation, but one cannot rely on its quantitative 
details in the present case. 

After this prelude on the density profiles for the well-studied case of polymer 
brushes grafted to planar substrates we turn to our results on the profiles in cylin- 
drical pores, Figs. ITHHT?)! While for short chains (such as N = 8, 12, 16) and low 
grafting densities these density profiles are typical of mushroom behavior, for larger 
iV one quickly reaches a state where the whole cylindrical brush is more or less 
uniformly filled with monomers. In this limit the density of end monomers reaches 
a maximum near the tube axis, but in comparison with the brushes on flat sub- 
strates (Fig. IT2j) the distribution of the end monomers is more uniform. Thus the 
approximation of Sevick [20] assuming a delta-function distribution of chain ends 
(like in the Alexander [UJ-de Gennes [TU] model) obviously is not accurate, at least 
for the parameter range accessible in our simulations. At this point we recall that 
Sevick [20] predicts that the brush height h decreases with increasing inverse cylin- 
der diameter D monotonously, while Milner and Witten [53] on the basis of their 
SCF theory suggest that the brush height h increases with 1/D (for small enough 
1/.D). The interpretation for the latter result is that in the cylindrical pore the blob 
diameters get smaller with increasing distance z from the grafting surface when 1/D 
increases, since the volume available for the blobs decreases, and hence the tendency 
of the brush towards stretching is more pronounced in the outer region of a brush in 
a cylinder than for a brush at a flat substrate. Of course, when h and D/2 become 
comparable, the opposite tendency of compressing the brush due to the confinement 
will outweigh this increase of h that may be present for small 1/D. 

Again we emphasize that the grafted monomer (adjacent to the wall of the tube) 
is at a fixed distance Ar = 1 inside of the atoms forming the walls of the tube, which 
are at the radial distance r = R from the tube axis. The sharp spike at r = R — 1 
in Figs. IT3k . [T4*b . [T5k . ITfik hence is always a representation of the delta function 
corresponding to the monomer adjacent to the rigid wall, at this fixed distance 
Ar = 1 inside of it. However, other monomers of the chain can come closer to the 
wall, when the chain bends back to the grafting wall. Therefore one finds a nonzero 
0(r) in between r = R — 1 and r = R, as is clearly seen in Figs. [THk and Hlk. 

When the chain length N is small enough (and/or the tube diameter D is large 
enough), both <f>(r) and p(r) are essentially zero for r — > 0. The pore center then 
is still free of monomers, and Figs. - EU then contain also the full information 
nn the distribution of the monomers as a function of the radial distance frnm, the. 



grafting wall, r', we simply can transform from r to z via r' = D/2 — r. However, 
when either N gets larger or D gets smaller (or both), these distributions <p(r) and 
p(r) get nonzero for r — > 0. This means that then distances larger than D/2 become 
possible {cf. Fig. I2^a)}. Then each value of r actually contains two entries for r', 
namely 

r' = D/2-r ,forr < D/2, and r' = D/2 + r, for r' > D/2. (25) 

The choice r' = D/2 + r has to be used when the monomer lies in the upper 
hemicylinder with respect to the grafting site taken as the lowest point of the lower 
hemi cylinder. Thus the tube axis in the simulation is at r' = R = D/2, since we 
arbitrarily defined r = R — 1 to be the position of the first monomer and put the 
wall at the radial distance r = R from the tube axis. 

Because of this ambiguity, it is also of interest to analyze separately the dis- 
tribution functions <fi(r') and p(r'). Fig. El gives an example (more details will be 
given in a forthcoming publication). Of course, for the cases which <f)(r — > 0) — > 
and p(r — > 0) the curves plotted in Fig. ITTk are identical to those of Fig. [T31 only 
now left and right is interchanged and the origin of the coordinate system shifted. 
However, this is not so in the cases where <f)(r — > 0) and p(r —>■ 0) are nonzero: we 
now see that <j)(r') and p(r') look qualitatively very similar to corresponding curves 
for brushes at flat walls. In particular, there is no singular behavior at the center 
of the tube, at r' — R all curves are perfectly smooth. Also the distribution of 
chain ends p(r r ) is perfectly smooth at this distance r' = R and carries considerable 
weight for distances r' > R if N is sufficiently large (and/or D sufficiently small), 
as Fig. [THb shows. 

It is also of interest to examine the distribution of the chain end p(y) as a function 
of the lateral coordinate y, parallel to the tube axis fFig. I18|). One can see that the 
width of this distribution does not depend much on the grafting density, and is 
largest for the smallest grafting density. No evidence for bimodal distributions, 
that are expected for narrow tubes and long enough chains due to the cigar-like 
configurations is found for the parameter range investigated, however. 

The profiles 0(r') can be used to calculate the first moment (r'), from (r') one 
obtains the brush height h, making use of Eq. f)23j) . Fig. EH shows a plot of the 
reduced height, h/N, versus N for D = 42 and D = 62 and two different grafting 
densities. A first important conclusion which one may draw upon inspection of Fig. 
Uni suggests that the brush height grows with increasing curvature of the tube, in 
contrast to the basic assumption in Sevick's theory j2Hj- Evidently, when the tube 
diameter diminishes the monomer density grows, cf . the inset of Fig. where also 
the monomer density profile for a flat substrate, that is, zero curvature, is displayed. 
Clearly, as the density grows, the excluded volume interactions "push" the brush 
height toward the tube axis. As shown in Fig. EE this effect is seen irrespective 
of particular grafting density a. At higher a, however, one can readily see that 
beyond a certain chain length (N = 32 for D = 42, and N = 40 for D = 62) the 
ratio h/N reaches a plateau, as expected on the ground of Eq. (J2J) whereas for the 
lower grafting density, o = 0.03, the onset of this scaling behavior probably occurs 
at much larger chain lengths N. 

6 Conclusions 

The configurations of flexible polymers endgrafted at the inner walls of cylindrical 
tubes under good solvent conditions were studied by two complementary approaches, 
namely by a scaling theory and by Monte Carlo simulation of a bead-spring model. 



The scaling treatment can identify the asymptotic power laws of the chain linear 
dimensions R gx , R gy and R gz and the brush height h in terms of the parame- 
ters of the problem, namely the chain length N, the grafting density a, and the 
tube diameter D. These power laws, however, are valid only asymptotically in 
the limit iV — > oo, keeping the scaled grafting density a = <riV 6//5 and the scaled 
tube diameter D = DN~ 3 ^ 5 fixed. In this limit, also the scaled linear dimensions 
Rga — RgaN~ 3 ^ 5 (a = x, y, z) and h = hN~ 3 / 5 will not exhibit any explicit depen- 
dence on N any more, but will depend on the scaled variables a and D only (see 
Fig. HI for example). 

It should be noted, that even in this limit the scaling theory is grossly over- 
simplified, since the different regimes that have been identified (mushrooms, cigars; 
compressed cigars; overlapping cigars; compressed brush; normal brush) are not 
separated by any sharp boundaries, but rather the different power laws that charac- 
terize these regimes merge gradually by smooth crossovers. These crossover regions 
may extend over several decades of the appropriate control variable (such as a). 
Therefore the lines in the "state diagram" of the polymer brush (Fig. ^ should not 
be mixed up with lines appearing in a thermodynamic phase diagram: here the lines 
only mark the center of the crossover region (fig. EJ) , and all the kinks in a scaling 
diagram such as Fig. 0] are rounded off. However, we are not aware of any analytical 
method that could provide the functions that describe these smooth crossovers. 

Real polymer brushes confined in such a cylindrical pore geometry will hardly 
obey these rather stringent conditions under which the scaling considerations are 
strictly applicable, and the same caveat holds for computer simulations. Neverthe- 
less the scaling approach is useful as a first guidance towards this rather complex 
behavior. 

The Monte Carlo simulations yield not only the gyration radius components 
of the chains, but also the full density profiles of the monomers, the distribution 
function of the chain ends, etc. It has been found that the chain ends are distributed 
all over the cylindrical tube, if the monomer density inside the tube gets large. While 
for brushes on planar substrates chain end distributions still have a broad maximum 
located roughly near the brush height, in the cylindrical geometry these distributions 
are much broader. It is also remarkable that this distribution spans the full range 
of radial distances (0 < r' < D), rather than being restricted only to the region 
from the grafting point to the cylinder axis (Fig. |2J), < r' < D/2. Our findings 
imply that theories that approximate the chain end distribution by a delta function 
5(r' — h) cannot be trusted. Thus, we identify cases where the brush height in the 
cylinder increases when the cylinder radius is decreased, in contrast to predictions 
of such theories. This means that for smaller D we find that the curvature slightly 
enhances chain stretching in the direction normal to the grafting substrate, but this 
effect is rather small. 

Due to technical difficulties of preparing well-equilibrated configurations of poly- 
mer brushes in cylindrical pores, only rather short chains could be studied, which 
do not allow a decisive test of the scaling theory yet. Nevertheless the scaling theory 
does provide a useful framework to analyze and discuss the Monte Carlo results. 

In future work, we plan to consider transport through such polymer-coated 
nanopipes. We also hope that our study will stimulate experimental work on such 
cylindrical brushes. Also a detailed comparison to self-consistent field treatments 
would be valuable. 

Acknowledgement: Support from the Deutsche Forschungsgemeinschaft (DFG) 
under project No 436BUL 113/130 is gratefully acknowledged. 




Figure 1: Blob picture of a polymer brush grafted on a flat substrate in the good 
solvent regime, (a) Alexander-de Gennes model: the chains are treated as linear 
"cigars" of blobs with uniform diameter £ = a^ 1 ^ 2 , a being the grafting density 
|3 EH]- (b) Non- uniform blob picture jH]: consistent with a decrease of the monomer 
density (j)(z) in the brush with increasing distance z from the grafting substrate 
surface, the blob diameter (which can also be interpreted as the screening length £ 
of excluded volume interactions, see (c)) increases with the distance z. 




Figure 2: (a) Choice of coordinates in a cylindrical tube of diameter D used for 
the analysis of the gyration radius components of the chains grafted on the interior 
pore surface. For each polymer chain a separate coordinate system is chosen with 
the origin at the grafting site whereby the z-axis is oriented perpendicular to the 
grafting surface and points towards the tube axis. The y-axis is parallel to the tube 
axis, and the x-axis is normal to both other axes in the tangential plane at the 
grafting site. Note that z-coordinates larger than D/2 are possible. This coordinate 
system is also used to study the distribution 0(f) of the monomers and p(f) of the 
free chain ends, (b) Schematic blob picture of a "cigar" , describing the state of a 
chain in a narrow tube in the dilute limit. 




Figure 3: Schematic diagram for the dimensions of a chain of length N in terms 
of tube diameter D and grafting density a. Straight lines indicate the crossovers 
between the various scaling regimes. For wide tubes (D > N 3 ^ 5 ) with increasing 
a only two crossovers are encountered: at a = a* = jV~ 6//5 from mushrooms to 
a quasi-flat brush, and then at a** = (D/N) 3 to a compressed brush. The latter 
regime ends at o = (1/4)D/N, where the tube is densely filled with monomers 
(<f) = 1). Since the line, corresponding to <ft = I (the regime > 1 is unphysical) 
meets the line a** at D = N, for D > N no compressed brush exists any longer. 
For narrow tubes, D < N 3 ^ 5 , instead of mushrooms one has cigars, elongated along 
the tube axis. Three sub-regimes need to be distinguished: For moderately narrow 
tubes, N 1 / 2 < D < iV 3 / 5 , three successive crossovers are distinguished: At a = a* = 
(ND 1 / 3 )^ 1 a crossover from swollen cigars to compressed non-overlapping cigars 
occurs. At a*** = _D~ 9 / 7 jV 3//r these cigars start to overlap strongly whereby for 
a *** < (j < a **** = jv 3 /} -7 their linear dimension in y-direction still exceeds the 
tube diameter. At or**** the chains behave similar as in a bulk semidilute solution of 
the same density and R gy = D. For o > a****, the regime of the compressed brush 
(as described for D > N 3 ^ 5 ) is entered again. For narrow tubes, iV 1/4 < a < iV 1/2 , 
the latter regime no longer exists, however. Finally, for ultranarrow tubes (1 < D < 
N 1 ^) only a single crossover from swollen to laterally compressed cigars is possible. 
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Figure 4: Schematic variation of the scaled linear dimensions R gy = R gy /N 3 ^ 5 , 
R gz = Rg Z /N 3 / 5 with the scaled grafting density a = criV 6 / 5 on a logdog plot. For 
wide tubes (D > 1, upper part) and moderately narrow tubes (D < 1 but 7J > iV 1 / 2 , 
lower part) there is no longer any dependence on chain length N. Since the crossovers 
at the various grafting densities a*, and a**** are not sharp, the dashed 

vertical straight lines should be understood as centers of crossover regions only, cf. 
text. 



Figure 5: Snapshot configurations of the brush in a tube of diameter D = 30 and 
length L = nD, for three choices of N: N = 8 (upper part), N = 20 (middle 

na.rtl and N = 9A flnwpr nartV Rnth virwr nf thp nntsirlp and nf rrnss spptinns ars 



10' 



DC io' 
V 



o N=16 
□ N=24 
o N=32 
N=40 
* N=52 
x N=64 




o o 



o o 



o o 



DC 

V 



o N=16 
□ N=24 




10 J 



1(T 



10" 



OC 

V 



o 0=0.0115 
■ o=0.1039 
1.43+/- 0.01 




10 



100 



Figure 6: Log-log plot of (R 2 Z ) (a) and (R 2 gy ) (b) versus grafting density a, using a 
perfectly flat substrate. The total number of chains ranges from 36 (cr = 0.011547) 
to 324 (a = 0.103926). Chain lengths N = 16,24,32,40,52 and 64 are included in 
the plot, as indicated. Straight lines indicate the theoretical power laws (R 2 gz ) oc a 2 / 3 
and (Rgy) oc cr -1 / 6 that apply in the semidilute regime of strongly stretched brushes 
under good solvent conditions (Eqs. (J2J), ©)• Upper part of (b) shows a plot of 
(Rgy) vs N, while (c) shows a plot of (R 2 ) vs. N. 
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shown as functions of the scaled grafting density criV 6 / 5 , in a log- log representation. 
Same data as shown in Fig. El are used, and again the theoretical power laws R 2 gz oc 
<t 2 / 3 and R 2 oc er -1 / 6 are indicated by straight lines. 
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Figure 8: Log-log plot of (R 2 gz ) (a) and (R 2 gy ) (b) versus grafting density a, for a brush 
in a tube of diameter D = 30. The total number of chains ranged from 144 (a = 
0.02) to 576 (a = 0.08). For a < 0.0356 chain lengths N = 12,16,20,24,32,40 
and 56 are included, while for larger a some of the larger chain lengths needed 
to be omitted since melt densities were reached inside the pore (recall Fig. |H1 the 
unphysical regime cannot be entered). 
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Figure 9: Same as Fig. |H1 but for D = 42. The total number of chains ranges from 
294 (a = 0.031) to 864 (a = 0.09). Chain lengths up to iV = 80 are included for the 
smallest and up to iV = 32 for the largest grafting density. 
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Figure 10: Same as Fig. |Hl but for D = 62. The total number of chains ranges from 
280 (a = 0.022) to 1144 (a = 0.09). Part (c) shows a log-log plot of (R 2 gz ) vs. N 
for two grafting densities. In parts (a), (b) we have included corresponding data 
for the flat brush on (R 2 gz ) (empty symbols) to show that at least for N = 32 and 
N = 40 the chains for the brush in the tube are slightly more stretched than for the 
flat brush. 
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Figure 11: Scaling representation where (upper part) and R 2 gy (lower part) is 
plotted versus a in double- logarithmic form, including data for various N and D. 
Open symbols: (thick small) D = 30, (thin big) D = 62. Full symbols: D = 42. 
Predicted power laws for the brush regime and the regime of the compressed brush 
are shown by straight lines. 
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Figure 12: Monomer density profile 4>(r') of a polymer brush on a flat substrate 
(a) and plotted versus the rescaled distance r'/N (b). All data refer to a grafting 
density a = 0.1 (N c h = 324 chains grafted to an area A = 3117.6). Seven chain 
lengths N are included in the figure, as indicated. The insert in the lower part shows 
a log-log plot of the brush height h versus chain length N. The distribution of free 
chain ends is shown in (c). 
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Figure 13: Radial density distribution function <f>(r) of the monomers as a function 
of the distance r from the cylinder axis, for a tube of radius R = 15, at a grafting 
density a = 0.02 (N c h = 144 chains), upper part, and the end monomer density 
distribution, lower part. Remember that at r = R — 1 = 14 the first monomer 
of each grafted chain is located, while the atoms forming the wall are located at 
r = R=lb. 
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;ure 15: Same as Fig. EH but for R = 31, a = 0.022 (N ch 




Figure 17: (a) Density profile <fi{r') of the monomers in the brush for D = 28 at the 
grafting density a = 0.02 and various chain lengths. Arrows point to the positions 
of the cylinder axis and the maximum distance from the grafting site that is possible 
in ^-direction inside the tube are marked by arrows, (b) Distribution p(r') of the 
chain ends, for the same parameters as part (a). 
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Figure 18: Probability distribution p(y) of the chain ends plotted vs. the in- 
coordinate along the tube axis, for D = 30 and two choices of N, N = 24 (a) 
and N = 40 (b). Several values of the grafting density are included, as indicated. 




Figure 19: Plot of the reduced brush height h/N in tubes with diameters D = 42 
and 62 as function of N at two different grafting densities. The inset shows the 
radial density distribution 4>{r') for iV = 32 and a = 0.03 for a flat substrate (the 
lowest 4>(r')), for a tube with D = 62 (intermediate), and for a tube with D = 42 
(largest 4>{r')). 
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